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ABSTRACT 

Collisions between giant molecular clouds are a potential mechanism for triggering 
the formation of massive stars, or even super star clusters. The trouble is identi¬ 
fying this process observationally and distinguishing it from other mechanisms. We 
produce synthetic position-velocity diagrams from models of: cloud-cloud collisions, 
non-interacting clouds along the line of sight, clouds with internal radiative feedback 
and a more complex cloud evolving in a galactic disc, to try and identify unique signa¬ 
tures of collision. We find that a broad bridge feature connecting two intensity peaks, 
spatially correlated but separated in velocity, is a signature of a high velocity cloud- 
cloud collision. We show that the broad bridge feature is resilient to the effects of 
radiative feedback, at least to around 2.5 Myr after the formation of the first massive 
(ionising) star. However for a head on lOkm/s collision we find that this will only 
be observable from 20-30 per cent of viewing angles. Such broad-bridge features have 
been identified towards M20, a very young region of massive star f ormation that was 
concluded to be a site of cloud-cloud collision by To rii et ah] (2011), and also towards 
star formation in the outer Milky Way by Izumi et al. ( |2014 ). 
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1 INTRODUCTION 


Understanding the galactic environment is an essential pre¬ 
requisite to understanding the Galaxy’s star formation rate 
(Fuji moto et al.||2014 Hughes et al.||2013 ). Interactions be¬ 


tween star-forming gas clouds are a major factor in cloud 
evolution, with events varying from tidal encounters to head- 
on collisions. While many of these are expected to be minor, 
leading to aggregation of material (Dobbs et al. 2015), more 


major collisions are being looked to as a possible formation 
mechanisms for super star clusters (SSCs) and massive star 


formation ( 

Furukawa et al.|2009a| Ohama et al.|2010 Fukui 

et al. 

2014 

). During such an encounter, shock compression 


at the collision interface leads to the rapid generation of 
unusually high densities compared to that found in an iso¬ 
lated, turbulent star-forming cloud. This high density slab 
can form massive cores either via an elevated Jeans mass or 
subsequent high accretion rates. Such a scenario tackles the 
traditional problem with forming massive stars, from which 
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strong radiation output during formation may terminate ac¬ 
cretion before a high mass is reached (McKee & Tan 2002). 
However, confirming these theories has been difficult due to 
the challenges of identifying a collision event. Specifically: 


(i) The collision itself will happen on relatively short 
timescales of typically a few Myr (e.g. it takes ~2Myr to 
cross a 7pc cloud at lOkm/s in the absence of braking). 

(ii) The number of massive star-forming regions in the 
Milky Way is low, implying that the collision frequency is 
similarly small. Observational and theoretical considerations 
suggest collisions between giant molecular clouds (GMCs) 
may happen once every 10 Myr, but the fraction of these that 
form stars will likely depend on morphology and collision 


speed (Tasker & Tan|2009 Takahira et al.|2014 Fukui et al. 
^0l5l. 


(iii) Massive star formation will result in feedback that 
can disrupt the cloud and the signature of a collision. Even 
if the signature is only partially disrupted, the system may 
just be studied as a massive star forming region, overlooking 
the possibility of a collision event. 
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(iv) The complexity of the multiphase interstellar 
medium (ISM) makes it difficult to differentiate between 
star-forming collision events and less extreme interactions. 
The situation is made worse given that most collisions are 
expected to occur in the densely populated spiral arms of 
the galactic disc ( Dobbs et al.||2013 Fujimoto et al.||2014 ). 


Recent work by e.g. Furukawa et al. (2009b); Torii et al. 
(2011) and Fukui et al.| (2014) has highlighted the potential 
importance of point (iii). They identify two distinct cloud 
components at different velocities surrounding SSCs. While 
relative velocities of the two clouds remove the possibility 
of a single bound system, their thermal properties suggest 
they are interacting with the star clusters. This raises the 
possibility that a collision between subregions of the two 
clouds may have been the trigger for the formation of the 
SSC. The SSCs themselves are observed to reside at the 
junction between the clouds, further supporting this idea. 


Recent numerical simulations by Dobbs et al. (2015) support 


the idea that collisions do occur over only small subsets of 
the clouds, leaving the extended gas to retain the kinematic 
structure of the pre-collision cloud, in qualitative agreement 
with the aforementioned observations. 

Collisions between GMCs have been studied using hy- 
drodynamical models both in smaller-scale dedicated colli¬ 
sion simulations (e.g. Inoue & Fukui 2013; Takahira et al. 


2014) and in the context of the entire galaxy (e.g. Dobbs & 


Bonnell|2006| [Dobbs et al.|2011||Tasker &; Tan|2009| |Taske7 


201 1[ |Dobbs et al.||2015| Fujimoto et al.||2014| . While only 

the former are capable of achieving the resolution required 
to follow the collapse into cores, the global simulations pro¬ 
vide a clue to the rate at which such interactions occur. 
They agree that interactions for a cloud should occur mul¬ 
tiple times per orbital period, a rate that could potentially 
drive the entire star formation rate of the galaxy if all colli¬ 
sions were productive ( Tan|2000 Fujimoto et al.|20l4 ). The 
local-scale simulations and observations, however, suggest 
collision results are velocity dependent. 

Recent N-body simulations by |Banerjee fe Kroupa| 
(2014 2015) found that the massive star cluster NGC 3603 


favours a monolithic formation mechanism over hierarchical 
assembly, the former of which is consistent with widespread 
massive star formation triggered by collisions (although not 
exclusively). 

Synthetic observations based on simple cloud collision 
models have also been tried. Keto & Lattanzio (1989) mod¬ 


elled collisions between uniform density high latitude clouds 
and then post-processed their results to produce optically 
thin 13 CO emission maps and line profiles. They found these 
could explain many of the observed features in high latitude 
clouds, including secondary peaks and broad wings. More 
recently, Duarte-Cabral et al. (2011) created surface den¬ 


sity position-velocity diagrams from hydrodynamical mod¬ 
els of colliding cylinders (which, seeded with a turbulent 
field quickly evolved into a collection of filamentary struc¬ 
tures). They found that their results provide a possible ex¬ 
planation for the trigger of recent star formation in Serpens. 
The clouds in this latter work were smaller than GMCs (the 
cylinders had initial radii of 0.25 pc and a typical length of 
lpc). 

In this paper we post-process a set of hydrodynamic 
and radiation hydrodynamic simulations, covering a range 


of possible scenarios including: non-interacting clouds coin¬ 
cident along the line of sight, a cloud that has evolved in a 
galactic environment, cloud-cloud collisions of different ve¬ 
locities and clouds with internal radiative feedback. Using 
these models we produce synthetic 12 CO position-velocity 
(p-v) diagrams which we use search for any signatures that 
are characteristic to a major cloud-cloud collision of the 
kind that might give rise to massive star formation. If one 
can determine the expected observational characteristics of 
enough theoretical models, any non-degenerate signatures 
should be very useful for interpreting real observations. 


2 NUMERICAL METHODS 

2.1 Hydrodynamic models 

2.1.1 Hydrodynamics simulations with ENZO 

Two of our hydrodynamic models of cloud-cloud collisions 


come from runs performed in 

Takahira et al. 

(2014 

), us- 

ing the grid-based hydrodynamics code, ENZO | 

Bryan et al. 


2014). The runs consist of two non-identical, turbulent, ide¬ 


alised clouds in a head-on collision. In one simulation the 
collision velocity is 3km/s and in the other it is lOkm/s. 
The large and small clouds involved in the collision have 
radii of 7.2 and 3.5 pc, initial turbulent velocity dispersions 
of 1.71 and 1.25 km/s and mean number densities of 25.3 
and 47.4cm -3 , respectively (i.e. the clouds have masses 417 
and 1635 Mq respectively). The clouds are initially spheri¬ 
cal and seeded with a turbulent velocity field. The collision 
itself is head on, in a frame such that the larger cloud has no 
bulk velocity with the smaller cloud moving towards it. We 
consider a snapshot of each simulation at the point where 
the maximum number of star-forming gas cores have formed, 
with a core being defined as gas within a contour of density 
of 0.3 x 10 4 cm -3 . The surface density of the 3 and lOkm/s 
collision snapshots is given in the middle and lower panels 
of Figure ^ respectively. Additionally, we consider one ex¬ 
tra snapshot from the lOkm/s collision model just prior to 
cloud contact, to isolate the signature of two clouds coinci¬ 
dent along the line of sight, but not interacting. We include 
this pre-collision snapshot in the upper panel of Figure [I] 
The limiting resolution (smallest cell size) in these simula¬ 
tions was 0.06 pc. 

A third hydrodynamic model involved in this compar¬ 
ison was of a star-forming cloud formed in a galaxy-scale 


simulation from Benincasa et al. (2013). The cloud was ex¬ 


tracted from the galaxy disc and evolved at higher resolution 
(Shima, Tasker and Habe, in prep.) to give a smallest cell 
size of 0.1 pc. The surface density distribution of this cloud 
is shown in Figure [2] 


2.1.2 Radiation hydrodynamic feedback models with 
SPH-NG 

To distinguish between observational signatures due to col¬ 
lisions and those due to stellar radiative feedback we also 
use 3 of the radiation hydrodynamic simulations of massive 
star feedback discussed in a series of papers by Dale et al., 


e.g. Dale et al. (2012) and Dale et al. (2013), namely models 


J, UP and UQ. These are smoothed particle hydrodynamics 
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Figure 1. Surface density plots of the ENZO cloud-cloud collision 
simulation snapshots. The upper panel shows the clouds just prior 
to a lOkm/s collision. The middle and lower panels show the 10 
and 3km/s collision models respectively. All snapshots are at a 
viewing angle perpendicular to the collision axis. 



Figure 2. A surface density snapshot of the cloud extracted from 
a galactic scale enzo simulation. 


simulations of radiative feedback from massive stars in tur¬ 
bulent clouds. They followed the gravitational collapse of 
a turbulent cloud and incorporated ionisation effects from 
any stars (sink particles) that attain mass ^ 20 M©. The 
effect of the ionising radiation is to heat the gas around the 
stars, resulting in hot ionised bubbles partially surrounded 
by molecular gas that is accelerated away from the ionising 
sources. Surface density plots of the snapshots we use are 
given in Figure [ 3 ] In model J the initial cloud is gravita¬ 
tionally bound, whereas in UP and UQ they are initially 
partially gravitationally unbound. 


2.2 Radiative transfer 


We use the TORUS radiation transport and hydrodynamics 


code to produce synthetic observations in this paper (Harries 
2000 Rundle et al.|2010| Haworth &; Harries|2012 ). Specifi¬ 


cally we use the non-LTE molecular line transfer developed 
by and discussed in great detail in Rundle et al. (2010) and 


summarised in Haworth et al. (2013). Given that it is docu 


mented comprehensively in these papers, we do not discuss 
the algorithm in any detail here. 

We map the results (densities, temperatures, velocities) 
of the various simulations onto the TORUS grid. The SPH 
data is mapped onto the TORUS grid using the method de¬ 
scribed in Acreman et al. (2010). We assume a fixed molec¬ 
ular abundance (8 x 10 -5 relative to hydrogen for 12 CO) 
and solve for the molecular level populations using a non- 
LTE statistical equilibrium calculation. We do include dust 
in these calculations and consider a standard ISM power 
law size distribution ( Mathis et al.||1977 ) with optical coef¬ 
ficients from Draine & Lee (1984) and a dust to gas ratio 
of 10 -2 . Once the level populations are converged they are 
used in further ray tracing calculations to produce synthetic 
position-position-velocity data cubes for as many observer 
viewing angles as desired ( Rundle et al.|[20 10). We use the 
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Figure 3. Surface density plots of the SPH-NG simulation snap¬ 
shots used in this paper. From top to bottom the panels are mod¬ 
els J, UP and UQ from Dale et al. White dots are sink particles. 



Figure 4. A schematic of the viewing angle convention used in 
this paper. A viewing angle of 6 = 7r/2, </> = 0 is along the collision 
axis. Note that the collision simulations take place in a frame in 
which the larger cloud is stationary. 

STARLINK software GAIA to collapse these data cubes along 
one spatial axis to transform them into p-v diagrams. 

We do not convolve the resulting p-v diagram to a beam 
representative of any real instrument since we are not trying 
to replicate a specific set of observations in this paper. The 
spectral and spatial resolution in the cloud-cloud collision 
models are 0.04 km/s and 0.2 pc. 

2.3 Summary of numerical models 

We post process a set of 7 numerical models to generate syn¬ 
thetic 12 CO J=l-0 p-v diagrams for clouds coincident along 
the line of sight, clouds undergoing collision, clouds subject 
to internal radiative feedback and a cloud that has evolved 
in a galaxy-like environment. For ease of reference and nav¬ 
igation of the paper we summarise the models in Table [l] 
Note that in section [6] we will also briefly consider the case 
of a collision between clouds, including internal radiative 
feedback. 


3 RESULTS OF SYNTHETIC OBSERVATIONS 

3.1 p—v diagram morphology 

We calculated 12 CO J=l-0 p-v diagrams for all of the nu¬ 
merical models summarised in Table[l] They are constructed 
by integrating along the entire x-axis at each velocity as a 
function of the y —axis on the image grid. The results are 
presented in Figures 5 through 8. We now discuss the differ¬ 
ent p-v diagram morphologies. 

3.1.1 Collisional models 

The p-v diagrams resulting from the collisional models 
(models A-C in Table [l]) are given in Figure [5] The up¬ 
per panel shows the p-v diagram for the two clouds in the 
10 km/s collision model at a time just prior to the colli¬ 
sion (model A) viewed along the axis of collision such that 
the clouds are coincident along the line of sight (0=7 t/2, 
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Table 1. Summary of the numerical models post processed in this paper. 


Model 

ID 

Scenario 

Origin 

Surface densiity 
snapshot 

p—v diagram 

A 

Turbulent clouds coincident along the line of sight 

Takahira et al (2014) 

Fig 1, top panel 

Fig 5, top panel 

B 

Turbulent clouds collided at 10 km/s 

Takahira et al (2014) 

Fig 1, middle panel 

Fig 5, middle panel 

C 

Turbulent clouds collided at 3 km/s 

Takahira et al (2014) 

Fig 1, bottom panel 

Fig 5, botom panel 

D 

A cloud that has evolved in a galactic disc 

Benincasa et al. (2013) 
(Shima et al. in prep.) 

Fig 2 

Fig 8 

J 

A gravitationally bound cloud with 
internal radiative feedback 

Dale et al (2012, 2013) 

Fig 3, top panel 

Fig 6, top panel 

UP 

A partially gravitationally unbound cloud with 
internal radiative feedback 

Dale et al (2012, 2013) 

Fig 3, middle panel 

Fig 6, middle panel 

UQ 

A partially gravitationally unbound cloud with 
internal radiative feedback 

Dale et al (2012, 2013) 

Fig 3, bottom panel 

Fig 6, bottom panel 


4 >=0 as defined in Figure Both clouds are clearly sepa¬ 
rated in velocity space by the known pre-collision velocity of 
10 km/s. Each individual cloud has a width in velocity space 
determined by their turbulent velocity dispersion. The low 
intensity, intermediate velocity gas is from material at the 
interface between the small and large clouds that has be¬ 
gun to collide/brake (the clouds are just touching at this 
snapshot in time). 

The middle and lower panels of Figure [5] show the 10 
and 3 km/s collision p-v diagrams respectively, also viewed 
along the collision axis (0— 7t/ 2, 0=0). Compared to the 
merely coincident clouds these diagrams are rather different. 
In the 10 km/s collision case, the two peak features are sepa¬ 
rated by lower velocity due to braking and are connected by 
a broad bridge feature. This signature is composed of the 
remnants of the two clouds (giving rise to the two peaks) 
and intermediate velocity gas at the interface between the 
two clouds (giving rise to the bridge). 

In the 3 km/s case two peaks are not discernible and 
there is therefore not any broad bridge feature. This is be¬ 
cause the clouds initially have a smaller velocity difference 
compared to the turbulent velocity dispersion of the clouds, 
making the bridge harder to distinguish, and this difference 
is rapidly reduced during the collision due to braking. The 
only signature different from that of an isolated turbulent 
cloud is the spike in negative velocity at ~ 0 pc, which comes 
from shocked gas. If the local turbulent velocity is lower it is 
possible that lower velocity collisions might be identified, for 
example Duarte-Cabral et al. (2011) have identified double 
velocity components separated by bridges in simulations of 
elongated collisions between low mass clouds at only 2 km/s. 

It therefore seems that a broad bridge separating two 
peaks and evidence of shocked gas are signatures of cloud- 
cloud collision, though whether this is identified is sensitive 
to the collision velocity relative to the turbulent velocity of 
the colliding gas. We will explore this velocity sensitivity 
further in section [5] 


3.1.2 Radiative feedback models 

The p-v diagrams resulting from the radiative feedback 
models (models J, UP and UQ in Table [l]) are given in Fig¬ 
ure [6] The panels are model J, UP and UQ from top to 
bottom. 

In the classical picture of an expanding Hu region, 


where there is a spherical shell of dense gas bounding the 
ionised gas, one would expect a circular or elliptical signa¬ 
ture in the p-v diagram. However in our results each di¬ 
agram consists of a single broad feature in velocity space, 
with the partially gravitationally unbound clouds consisting 
of multiple components along the spatial axis. This signa¬ 
ture arises rather than the elliptical one because these Hu 
regions are actually very leaky, losing up to 95 per cent of 
the ionising photons (this has also been suggested observa- 
tionally, by e.g. Beaumont & Williams][ 2010). Only a small 
fraction of 47r centred on the ionising sources is subtended 
by dense molecular material and so only a small fraction of 
the classical sphere is actually accelerated. 

The p-v diagrams also show small spatial scale, high 
velocity features, which are globules (interior to the HII re¬ 
gion) or clumps of gas at the edge of the HII region accel¬ 
erated by the ionising radiation field. We illustrate this in 
Figure [ 7 ] by showing the integrated emission of model UP 
with contours from the highest velocity channels overlaid. 

Overall these p-v diagrams of radiative feedback mod¬ 
els show an imprint of the initial turbulent velocity structure 
coupled with relatively localised components of the cloud ac¬ 
celerated by the ionising radiation field. The p-v diagrams 
for models with radiative feedback unsurprisingly show sig¬ 
natures of high velocity gas, however none show the broad 
bridge feature that we find in our 10 km/s model, supporting 
its use as a signature of cloud-cloud collision. 


3.1.3 A larger cloud from a galactic disc 

Figure |8] shows the p-v diagram for the more complicated 
cloud taken from galaxy scale simulations of[Benincasa et al.| 
( |2013| ). Of course a single cloud is not really representa¬ 
tive of the whole range of clouds evolving in a galactic disc 
(probing this full range is beyond the scope of this paper), 
however we include one to compare with a more complicated 
cloud structure. This cloud has evolved in a Milky Way-type 
galaxy without a grand design spiral where the environment 
is highly dynamic, it is therefore complex in geometry due 
to gravitationally driven mergers and tidal events happen¬ 
ing regularly. This cloud is also much larger than the other 
clouds considered here. 

The p-v diagram is highly disordered and very broad 
along both spatial and velocity dimensions. The different 
features predominantly come from the series of filamentary 
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Figure 5. 12 CO J=l-0 position-velocity diagrams of the ENZO Figure 6. 12 CO J=l-0 position-velocity diagrams for the feed- 

models of cloud-cloud collisions. The upper panel is for a sim- back simulations of Dale et al. From top to bottom the panels are 

ulation snapshot just prior to the collision (the clouds are just models J, UP and UQ. 

touching) in the lOkm/s collision model, with both clouds along 

the line of sight. The middle and lower panels are from the 10 and 

5 km/s collision models respectively at the point of maximum core 

formation in each model (c.f. |2.1.l] >. 
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Figure 7. Integrated CO J=l-0 emission from model UP with 
contours overlaid from the highest and lowest velocity channels. 
High velocity features are globules within the HII region and ac¬ 
celerated gas at the H il region boundary. 


structures which are attached to the main high density cen¬ 
tral clump. At ~-50pc there is a feature that might even 
be interpreted as a broad bridge, however this is quite mis¬ 
leading given the scale of the model, as there is actually a 
spatial offset of ~ 10 pc between the two intensity peaks. Fur¬ 
ther inspection of the full data cube (and Figure [l]) reveals 
this these separated velocity features are due to multiple 
filaments almost coincident along the line of sight. 

This system highlights that p-v diagrams require coor¬ 
dination with other techniques to robustly interpret observa¬ 
tions. It also demonstrates that at larger spatial scales inter¬ 
pretation will be more difficult since there will be a greater 
combination of processes contributing to the emission. Ide¬ 
ally one wants to observe on the scale of the bulk dynamics 
of one specific process (such as feedback or a collision) with 
only underlying turbulence as an additional factor, in addi¬ 
tion to the larger scale observations. 



-100 -50 0 50 100 

Position, pc 



Position, pc 

Figure 8. 12 CO J=l-0 position-velocity diagrams for the com¬ 
plex GMC from a galactic scale model. The upper panel shows 
the whole p-v diagram and the lower zooms in on the region at 
~-50pc to illustrate that there is no broad bridge feature since 
the intensity peaks are spatially offset. 


4 COMPARISON WITH OBSERVATIONS OF 
M20 

We find that the broad bridge feature is a signature of cloud- 
cloud collision that does not arise due to radiative feedback 
or for clouds merely coincident along the line of sight. Iden¬ 
tifying such a feature towards regions of star formation (in 
particular massive star formation where radiative feedback is 
prevalent) may therefore provide evidence of star formation 
triggered by cloud-cloud collision. In this section we will 
discuss 12 COJ=1-0 observations towards M20 taken with 
Mopra, which show candidate broad bridge features. 


4.1 Details of the observations 

Our 12 CO J=l-0 observations towards M20 were carried 
out using the 22m ATNF Mopra telescope in Australia in 
2011 October. The Mopra backend system “MOPS” pro¬ 
vided 4096 channels across 137.5 MHz in each of the two or¬ 
thogonal polarizations, and the corresponding velocity reso¬ 


lution and velocity coverage were 0.088 km/s and 360 km/s, 
respectively, at the frequency of 12 CO J=l-0 115 GHz. The 
OTF (on-the-fly) mode was used, and the pointing accuracy 
was checked every OTF scan to be better than T^with SiO 
maser observations at 86 GHz. The typical system temper¬ 
ature was about 500 K at 115 GHz. The absolute intensity 
calibration were made with the observations of Orion-KL 
(R.A., Dec.) = (5:35:14.5, -5:22:29.6) by comparing the re¬ 
sults of Ladd et ah] ( |2005 ). The obtained spectra were grid- 
ded to a 15"spacing and then were spatially smoothed to a 
beam size 45^. The achieved rms noise level is ~0.2K per 
channel at velocity resolution 0.9 km/s. 


4.2 Discussion of M20 observations 

Integrated 12 CO J=l-0 emission over two velocity ranges are 
overlaid upon an optical image in Figure [9] The black cross 
represents the exciting O star. In this paper we only discuss 
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Figure 9. The two colliding clouds identified by Torii et al. (2011) are presented with the Mopra 12 CO J=l-0 data at dh'Tesolution. The 
background is an optical image of M20 (Credit: Todd Boroson/NOAO/AURA/NSF). The cross indicates the existing O star, responsible 
for the ionising radiation field. Dashed lines show the regions used for the position-velocity diagrams in Figure [To| 


p-v diagrams over the white boxed regions, the full presen¬ 
tation of the Mopra dataset will be published elsewhere. 

shows our 12 CO J=l-0 position-velocity dia¬ 


Figure 10 


grams towards M20. The upper panel is for gas away from 
the exciting 07 star in M20 (the right hand white boxed 
region in Figure [9j whereas the lower panel is for material 
close to the star (the left hand white boxed region in Figure 
[ 9 ]). The 07 star is at a declination of 23h 031 x . Assuming a 
distance of 1.7kpc to M20, both bridge features have a size 
of roughly a few parsecs (similar to our simulation results). 


There is striking morphological similarity between the 
observations in Figure and our simulations in Figure [5] 
though in the lower panel of Figure [T0| the intermediate ve¬ 
locity gas looks like it is at least partially spatially offset. 
This could imply chance spatial correlation between three 
velocity components, or might be due to the action of the 
ionising radiation field from the O-star disrupting the colli¬ 
sion signature. 


M20 is a very young site of massive star formation that 
was concluded to be the site of a cloud-cloud collision by 


Torii et al. (2011). The numerical models in this paper cou¬ 


pled with these MOPRA observations support this hypothe¬ 
sis. Given that this system is very young, we speculate that 
either the signature of the collision has not yet been dis¬ 
rupted by feedback (though there could be signs of feedback 
in the broad bridge close to the O star), or that the colli¬ 
sion is ongoing between different components of the parent 
molecular clouds. 


A broad bridge feature was also recently identified in 
p-v diagrams towards star formation in the extreme outer 
galaxy by Izumi et al. (2014), which they propose was trig¬ 
gered by a high velocity collision. Our simulations support 
their conclusion. Galvan-Madrid et al. (2010) also identify 


broad bridges in the W33A high mass star forming region, 
which they interpret in terms of converging flows. 

In general broad bridge features can only really be used 
in conjunction with other diagnostics. For example, colli¬ 
sions between extended clouds should result in undisturbed 
subsets of the clouds that will be observable as spatially seg¬ 
regated red and blue shifted gas (e.g. Torii et al.|2011 Fukui 
et al.||2014 ) 


5 VIEWING ANGLE SENSITIVITY 


So far the discussion of p-v diagrams for collisional models is 
all regarding a viewing angle along the collision axis (6— 7r/2, 
(j )—0 in Figure [4]). 

For the simple geometry of a head on collision we can 
make a rough estimate of the fraction of viewing angles that 
a broad bridge is discernible. During the collision the clouds 
are travelling at velocities v\ and V 2 (which may be slower 
than the pre-collision velocity due to braking) with turbulent 
velocity dispersions Avi and Av 2 . If we require the clouds 
to be separated by more than their velocity dispersions then 
simple geometric consideration under the convention given 
by Figure [4] yields 

(Aiq + At> 2)/2 < sin(0) cos (4>)(vi — V 2 ) (1) 


as the criterion for identifying the broad bridge in a p-v 
diagram. Note that we have assumed that the velocity dis¬ 
persion is independent of viewing angle. Integrating over 0 
and (j) using Av\ — 1.71, Av 2 = 1.25 (c.f. section 2.1.1) and 
V 1 —V 2 ~ 3.5 km/s (by visual inspection of the middle panel 
of Figure [ 5 ]) gives an estimate of 28 per cent of viewing an¬ 
gles over which the 10 km/s collision (braked to 3.5 km/s) 
will be identified. 
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Figure 10. MOPRA 12 CO J=l-0 position-velocity diagram ob¬ 
servations of candidate broad bridge features towards M20, a 
young site of massive star formation. 


To check our simple estimate we produce p-v diagrams 
over a range of different viewing angles. To summarise the 
results of this process we calculate the mean profile over the 
middle of the p-v diagram where the bridge feature is found. 
For the lOkm/s collision model this is shown for a selection 
of example 0 in the top panel of Figure [TT) At some viewing 
angles two peaks in the p-v diagram are discernible and we 
conclude that the broad bridge can be identified, whereas at 
others there is a single (albeit broad) peak and the broad 
bridge is not identified. We create a grid of viewing angles 
spanning 6 = 0 to 90 in 10 degree intervals and <f> = 0 to 
180 in 20 degree intervals and assume that the symmetry of 
the problem allows us to translate to other viewing angles. 
We then assume that if there are two peaks in the averaged 
line profile then the broad bridge (and hence the collision) 
is identified. Using this definition and spanning our grid we 
estimate that only 20-30 per cent of viewing angles over 4n 
steradians will observe the broad bridge feature in the p- 
v diagram (at least at this collision velocity and geometry 
of a head on collision), consistent with our simple estimate. 



Velocity, km/s 

- 1 - 1 - 1 - 1 - 1 - 

clump 1 - 

clump 2 . 



Velocity, km/s 

Figure 11. The upper panel shows synthetic P-V diagram pro¬ 
files averaged at each velocity, illustrating the presence of two 
peaks at some, but not all, viewing angles (where viewing angles 
are specified in degrees using the convention of Figure [4] where 
4 > = 0 and 0 varies). Note that this is just an illustrative selection 
of viewing angles and the profile is reasonable symmetric when 
the viewing angle is translated by 180 degrees. The lower panel 
shows the same averaging over the two broad bridges observed 
towards M20 (clump 1 and 2 are the two regions in the upper 
and low panel of Fignre [l0| respectively). 

Observing a broad bridge feature might therefore provide 
evidence for a collision, but because they are only identi¬ 
fied from a small range of viewing angles not observing one 
is not strong evidence for a lack of collision. Conversely, 
should broad bridge features become commonly identified, 
this would imply either that collisions take place at veloc¬ 
ities much greater than the turbulent velocity (and so are 
observable at more viewing angles) or there is some further 
mechanism at work that also gives rise to this signature, 
making it degenerate. 

Note that these conclusions regarding the fraction of 
viewing angles from which the collision is visible are specific 
to the geometry of head on colliding clouds with no extended 
structure that does not undergo collision. It is also specific 
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Figure 12. A 12 CO J=l-0 position-velocity diagram of a new 
ENZO model of a cloud-cloud collision including star formation 
and radiative feedback, to be published in Shima et al (in prep). 
At the time in the simulation that this diagram is produced there 
are many ionised bubbles due to radiative feedback. This illus¬ 
trates that the broad bridge is resilient to the effects of radiative 
feedback. 


to the turbulent spectrum of the clouds. Collisions between 
clouds with lower turbulent velocities will be easier to iden- 
tify. 

The bottom panel of Figure EH shows the mean profile 
over the middle of the p-v diagram for the observations to¬ 
wards M20, clearly showing two discernible peaks similar to 
those that we see in the models. The M20 peaks are sepa¬ 
rated by a larger velocity than our simulations, which could 
mean that the collision was between higher velocity clouds 
or that it is at an earlier stage than our snapshot times and 
has not undergone much braking. If the latter, given that 
stars have already formed, it is possible that the collision is 
taking place sequentially along different components of the 
two clouds. 


6 MODEL LIMITATIONS 

The cloud-cloud collision models considered two isolated, 
initially spherical, turbulent clouds undergoing a head on 
collision. In reality these initial conditions are not neces¬ 
sarily representative of those formed in galaxy scale models 
(e.g. Rey-Raposo et al.|2015 ). It would be useful to perform 
a similar analysis using two clouds derived from galaxy scale 
models, and using different impact parameters. Even better 
would be to re-run collisions which happen “naturally” in 
galactic scale models at higher resolution. Ideally future sim¬ 
ulations would also model the formation of any star(s) using 
sink particles. This would also allow us to search for segre¬ 
gated red and blue-shifted cloud components about the star 
cluster, which are interpreted as evidence of cloud-cloud col¬ 


lision by Furukawa et al. (e.g. 2009b); Torii et al. (e.g. 2011); 


Fukui et al. (e.g. 2014). 


We have also not determined how long this signature 
survives after the collision. Given that internal radiative 
feedback seems to be primarily accelerating small clumps 
(see section 3.1.2), it is conceivable that the broad bridge 
might be long lived. To gain an initial insight into the effect 
of radiative feedback when it comes to disrupting the broad 
bridge we postprocess a snapshot from the ongoing simu¬ 
lations of Shima et al (in prep). These are similar to the 
cloud-cloud collision models in this paper, but the clouds 
are larger, colliding at 20km/s and the simulations also fol¬ 
low the formation of stars and include radiative feedback. 
For now, we study one snapshot 3Myr after the onset of 
collision and 2.5 Myr after the formation of the first massive 
star, at which time there are multiple ionised bubbles. We 
plot the p-v diagram of this snapshot in Figure [12] in which 
a broad bridge is definitely discernible. At least in this sce¬ 
nario, at this snapshot in time, the broad bridge is still ob¬ 
servable despite radiative feedback. We will study the time 
evolution of the broad bridge in simulations with radiative 
feedback in subsequent work. 

Since we are working on the premise that the cloud 
collision triggers the formation of massive stars, it should be 
unlikely that a site of cloud-cloud collision will be subject 
to an external ionising radiation field unless the collision is 
sequential along the cloud length (i.e. the collision is ongoing 
after the formation of massive stars) which is not the case 
in the collision models studied here. 

Finally, the astrophysical scenarios that we explore may 
not be exhaustive. For example, we do not investigate p- 
v diagrams for clouds disrupted by supernovae, large scale 
flows due to gravity, radiative feedback in less leaky HII re¬ 
gions (which should give the elliptical signature discussed in 
section 3.1.2) or those traversing spiral arm shocks. In par¬ 


ticular outflows can produce high velocity features, which 
we do not study here. If any of these processes can also give 
rise to a broad bridge then the signature becomes degener¬ 
ate (though potentially still useful in conjunction with other 
diagnostics). 


7 SUMMARY AND CONCLUSIONS 

We have calculated synthetic p-v diagrams for a number 
of different astrophysical systems: cloud-cloud collisions, 
non-interacting turbulent clouds coincident along the line 
of sight, turbulent clouds with internal radiative feedback 
and a GMC evolving in a galactic disc. We compare these 
to try and identify characteristic signatures of cloud-cloud 
collisions. We draw the following conclusions from this work: 


1) Cloud-cloud collision models give rise to a broad 
bridge structure in p-v diagrams. This broad bridge is 
two intensity spikes separated by lower intensity emission 
across the velocity axis. This feature is not reproduced by 
p-v diagrams from any of the other scenarios, potentially 
making it a useful signature for identifying cloud-cloud 
collisions. We also find instances of this broad bridge feature 
towards M20, a very young site of massive star formation 
which was concluded to be a site of cloud-cloud collision 
by Torii et al. (2011). Our models therefore support the 
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conclusion that M20 is a site of cloud-cloud collision. 

2) The broad bridge feature is observable in our lOkm/s 
collision model, but not in our 3km/s model. We conclude 
that in order to observe the broad bridge, the difference in 
line of sight velocity between the two colliding clouds at the 
time of observation has to be greater than half the sum of 
the turbulent velocity dispersions of the clouds. 

3) Using the criterion from conclusion 2, we estimate 
that for the specific case of our lOkm/s head on collision 
model, the broad bridge is only observable over 20-30 
per cent of 47r steradians. Given that it is also not clear 
how long the broad bridge survives before being dis¬ 
rupted by star formation and/or feedback, not observing 
a broad bridge can therefore not be used to rule out 
cloud-cloud collisions. Conversely widespread identification 
of this feature might suggest low turbulence, fast colli¬ 
sions, the importance of impact parameter, or alternatively 
that some other mechanism also gives rise to a broad bridge. 

4) Using ongoing models from Shima et al (in prep), we 
preliminarily show that the broad bridge feature is resilient 
to the effects of radiative feedback at least up to 2.5 Myr 
after the formation of the first massive star. 

5) In general, a broad bridge feature seems to provide a 
useful signature of cloud-cloud collision, but should be 
used in conjunction with other diagnostics such as channel 
maps and checking for high velocity segregated blue and 
red shifted clouds near the proposed collision site that may 
be relics of the pre-collision clouds. 
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